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I. INTRODUCTION 



The spin-isospin properties of nuclei are a long standing and still very interesting subject 
of nuclear physics New theoretical and experimental developments are seen in the 

study of nuclear spin-isospin response functions over the last decade. They promoted many 
interesting problems and active researches in this fields 0. 

In the analysis of the responses in the quasi-elastic region, nuclear currents with the spin 
and isospin degrees of freedom are usually separated into the isovector spin-longitudinal and 
-transverse components. They are characterized by the operators, t(<t • q) and t{(t x q), 
respectively, where q is the transferred momentum to the nucleus and q = q/q. The spin- 
transverse response function Rrp{q,uj) has long been observed in a wide range of q and the 
transferred energy uj by electron scattering However, hadronic probes are needed to 

study the spin-longitudinal response function c<j). 

Owing to the great progress of the experimental technique, complete measurement of the 
polarization transfers Dij was first carried out at LAMPF for the quasi-elastic {p,p') 
scattering, from which they extracted the ratio R^^/Rrp by use of an eikonal approximation. 
Measurement of {p, n) reaction is more difficult but more preferable because it can exclusively 
extract the isovector part. Recently LAMPF group |llO|-|T2| observed Dij of the quasi- 
elastic reaction and extracted R^/Rj. in a similar way. An interesting finding of 
these experiments is that the ratio is close to or less than unity in the quasi-elastic region. 

This was surprising because the standard nuclear model predicted that the ratio was 
much larger than unity at relatively low u. Based on the random-phase approximation 
(RPA) in nuclear matter, Alberico et al. pointed out that R^^ is enhanced and its peak 
position is shifted downwards (softening), whereas Rj. is quenched and its peak is shifted 
upwards (hardening) around the transferred momentum q = 1.75 fm~^. Therefore the ratio 
Ri^/Rrp extremely exceeds unity especially at relatively low uj. However, such behavior has 
hardly been seen in the experiments. 

This contradiction has been challenged from various aspects, such as finite size effect of 
nucleus jl^HlTi; relativistic RPA approach |jl8|, nuclear correlations beyond RPA |l9|-p3 



effects of absorptions and distortions ||T7|j2^ -p7| , etc.. Among them here we investigate effects 
of A-mixing and dependence on the nucleon particle-hole (ph) and A-hole (Ah) effective 
interactions. 



The analysis of Alberico et al. |]T3[ took account of the A degree of freedom and used the 
(tt + p + g') model for the effective interaction, namely one-pion exchange -|- one-rho-meson 
exchange + the contact interaction specified by the Landau-Migdal parameters g'^s. There 
appear three different ^'"s relevant to the interactions between ph and ph, between ph and 
Ah and between Ah and Ah. They are denoted by g'j^^^, g'^/^ and (^aa^ respectively. 

For computational simplicity, most of previous works adopted the universality ansatz 
o'nn — 9na — 9'aa- However it has no theoretical foundation nor was supported by 







various estimations 0. For instance, the phenomenological analysis yields g'j^^ ~ 

0.6 ~ 0.7. G-matrix calculations suggested g'^^ ~ 0.5 and g'j^^ ~ 0.4, but the induced 



interaction increases and the finite size effect decreases them [p2| -p^. A recent estimation 
by Brown et al. [^] gives g'^^^ ~ 1.0, g'j^^ ~ 0.33 and (^^^ ~ 0.5. 

Considering this situation, we developed a RPA formalism for the finite nucleus composed 
of and A, in which the universality ansatz is not adopted. Then we investigated the effects 
of A and the g' dependence of i?^ and Rrp. Similar analysis was performed in the Fermi gas 



model by Shiino et al. However the model gives unreasonably large enhancement and 



softening of i?^ and hence it is somewhat misleading. Here we perform a more realistic 
and detailed analysis for finite nuclei. 

The response functions including A^ and A constitute of three different components, 
]^[NN]^ j^[NA]^_ ^[AAfj-j g^^^ R^^^\ which Correspond to the processes depicted in Fig. |l|. 
j^i^N] represents the process in which the nuclear current first creates a ph state and anni- 
hilates it at the end. i?'^^] corresponds to the process which start with the Ah creation and 
end with the ph annihilation, and R^^'^^ corresponds to the inverse one. R^^^'^ represents 
those which start from the Ah excitation and end with its annihilation. We call them the 
process-decomposed response functions. 

We will see that the Ah configuration is crucial for the enhancement of R^ and plays 
some role for the quenching of Rj. for ordinary interactions. These effects come mainly from 
which is sensitive to g'jy^- This indicates that we must avoid the universality ansatz 
and treat g"s independently. 

In many of the previous works with the universality ansatz, the decomposition did not 
appear explicitly. The response functions have implicitly been calculated by 



fN 

where /a? and /a are the vrA^A^ and ttA^A coupling constants, respectively, and /a = S.O/at 
is usually used. However, observed response functions depend on probes. For instance, the 
isovector spin-transverse response function Rp^ ■* observed by (e, e') is approximated by 

p(e,e') _ r^[NN] ^ f^NA p[JVA] / f-yNA \ ^ p[AA] /„x 
J-yNN \ J^NN / 

where f^^j^f and f-yNA specify the isovector magnetic coupling strength of 7A^A^ and 7A"A 
vertices, respectively. We also check this coupling constant dependence. 

In Sec. H we present a formalism for calculating the response functions under the ring 
approximation in the finite system which is composed of A^ and A. The finiteness is handled 



by the continuum RPA with the orthogonality condition [40,17,41|. The universality ansatz 



for g' is removed. In Sec. |I| we show the effects of the A mixing on the response functions. 



and analyze them in the form of the process-decomposed components. The energy-weighted 



and non-weighted sums are also discussed. In Sec. |IV| we investigate the effective interaction 
dependence of i?^ ^ and R^i^ ■ We compare the results obtained by various values of g' . We 
also show some results with lighter p-meson effective mass m*. In Sec. |V] we present our 
calculation of R^^'^ with experimental data. Its dependence on g' and the ratio f\NA/ f:!/NN 
is shown. The summary is given in Sec. [V^. 



II. FORMALISM 

In this section, we formulate a method to calculate the response functions for a finite 
system composed of and A in the ring approximation. For simplicity, we only consider 
doubly (sub-)closed shell nuclei, so that the particle and the hole states are well defined and 
the spin of the ground-state is zero. 

We express the spin and the isospin operators of (cr and r) and the transition operators 
between and A {S and T) in unified spherical tensor forms, a^J-"-^^ and 

0-°^''^^ = 6ab for s = 0, 

- - S,, - i-rS-,^ for . = 1, 

with a, 6 = or A, and in a similar way for r**^"^). 



A. Spin-isospin polarization propagator 



We define the spin-isospin (transition) current operators as 



■tu{ab)( \ — 



' SfJ, 



k=l 



tu(ab) I 



r;rk 



(4) 
(5) 



where is the position vector of the fc-th particle in the nucleus and j*^"^-* operates only on 
the type (6) particle. They are separated into the angle and radial parts of r as 



iJM m 



JisJM V ) 



(6) 



where 



.tt/(a6)/ X 
JesJM V J 



A 
k=l 



jl:fSir;rk) ^ {rr^)^ [^'Y\n.J ® (a^^'^^O 



J 

M 



(7) 
(8) 



A 



We then introduce the spin-isospin polarization propagator as 
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where j = j — ( \[^o I j I ) is the current fluctuation and | ) and Eq denote the ground- 
state of the nuclei and its energy, respectively. The second identity comes from the assump- 
tion that the total angular momentum of the ground-state is zero. 

In the pure shell model without the residual interaction, the total Hamiltonian H is 
replaced by the uncorrelated one H^^\ the sum of the single-particle Hamiltonian. Then the 
uncorrelated polarization propagator is given by 



n<ir^(r,r» = 5.AeE 
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where | $o ) is the ground state of and \p{A)h) denotes a ph{Ah) state. Here e^^ 



denotes the single-hole energy and e^?]^-, is the single-particle energy of the type (a) particle. 



In a continuum RPA [Q, the sum over the discrete and continuum particle states is 
carried out by use of a single-particle Green's function, 

1 



(r 



e — + if] 



r') 



(12) 
(13) 



where h^""^ is the single-particle Hamiltonian of the type (a) particle and y^sjjyii^ 



X 



with the spin function Xm<> and the isospin function r/^t 



Using g^^\r, r'; e) and the single- hole wave function u^p{r)^ we get p^JT7[ 



<^ {n(3)eocc 

^aNB^,) ^(/5,a)— — g)^\r,r;u} + e^p) B^.Jj/ {/3,a) 
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+SbNB: 
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where a = (iaSajata^a) — i.^p^pjp'tiB''^^)^ and 



esJ ' 
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for s = 

for s = 1, Si = S2 = 1/2 

for s = 1, si = 3/2, S2 = 1/2 

for s = l, si = 1/2, S2 = 3/2 



(15) 



(16) 



B. Ring approximation 

We take into account the nuclear correlation by the ring approximation. We assume that 
the effective interaction is the charge-independent local two-body force, which is written as 

ss'W JMtu ^ 

(17) 



Then the polarization propagator satisfies the ring equation 



nlv.'(^,^» = n(°^,:,v(^,^» 



J2 I '"i^(^nr2^cir2 n(°)*i7£i.i('^,n;o;)iyy^^,^^2S2('"i.^2;w) 



xn 



.,,(r2,r», (18) 



where 11 and W are the matrices with respect to the particle types, 

j^{NN,NN) j^{NN,N/^) j^{NNAN) 
n = Jl{N/^,NN) -Q(ArA,AfA) -Q(iVA,AAr) 
-Q(AAr,7VAr) -Q(AAf,AfA) -Q(AAr,AAf) 

y^(NN,NN) y^{NN,N/:^) y^{NNAN) 
W = yf/(A^A,JViV) lYiNA,NA) p^(iVA,AJV) 
^r{AN,NN) p^(AAf,7VA) p^(AJV,AiV) 



(19) 



(20) 



Using the symmetry of ly, we introduce the grouped notation lyl"^' for W^"''''^'^^ as 

^[iViV] ^ y^^(NN,NN)^ (21) 

lylNA] ^ y^[AAf] ^ p^(7VAr,7VA) ^ y^(7VAr,AiV) ^ y^[NA,NN) ^ p^(A7V,Ar7V) ^ (22) 

p^[AA] ^ p^(AfA,iVA) ^ p^(7VA,AAf) ^ ^^{AN,NA) ^ p^{AiV,A7V)^ (23) 

and correspondingly the grouped representation of 11 as 

jjlNN] ^ Yl(^N,NN)^ (24) 

Yl[NA] ^ -QCAfAT.ATA) _^ ]-[(AfAr,AAr) (2^^ 

]-[[AJV] ^ -(-[(iVA.TVJV) _^ ]-[(AAf,7VJV)^ (2g^ 

■q[AA] _ jj{NA,NA) _|_ ^-[(AfA.AAr) _|_ ■Q{AAf,ArA) _|_ ^-[(AAr.AAf) (2j-j 

Then we can reduce the dimension of 11 and W in Eq. (|l5), and express them as ||39|| , 



n = 



w = 



lYlAN] ^,[AA] 



Yl[NN] YliNA] 
Yl[AN] jj[AA] 

The ring equation ( ]T8| ) is now symbolically written as 

n = n(°) + n(°)vrn, 

with 



We also get the symmetry, Ht^^] = U^^^\ since W^^^'^ = W^^^^ 



(28) 



n(o) 





n(o)[AA] 



(29) 
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C. Spin-longitudinal and -transverse modes 

In this paper we restrict ourselves to the response functions for the isovector spin- 
dependent currents {s = t = 1). So we use the abbreviations, 

■ i'(af)) ■t=li'{ab) -i-i-!^(afe,cd) -i-i-t=14'=l!/(afe,cd) /'Q-l\ 

JUM =3ls=lJMi ^^JW = ^^,}ts=U's'=l ■ \^^) 

We define the spin-longitudinal and -transverse currents, jl and jt-, in the momentum 
representation as 

j^^('^^)(q)^^j^^('^^)(q;r,), (32) 

k=l 

jTj:"'\q)^j:jTt'\q-^r,l (33) 
fc=i 
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where 
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They are separated into the angle and radial parts of q as [|14 



UM 



UM 



where 



FljMfii^q) = E 
KQ 



Alii ■ h 



JKl\ 



\M- 



I J 


K 


1 ^ 














-Q 


-/i j 





with 



aji = {J010\m), 



bjKe = jQi2J + l){2K + l) aKi 



1 1 1 

J K e 



(34) 
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(40) 
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The current j{q] r^) is related with j(r; r^) defined in Eq. (P) as 

feljuiQ^rk) = f r^dr j';^j'ij{r;rk)je{qr), (42) 
Jo 

with the £-th order spherical Bessel function ji{qr). 

In terms of these currents, we define the spin-longitudinal and -transverse polarization 
propagators as 



nu(ab,cd) I I \ 
T^^r '{q,q;uj) 



Jl iq) 
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~ v{ab) . , 
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1 



From Eqs. (|36|) , (|37D and (|^), they can be rewritten as 



jtu [q) 



(44) 



u(ab,cd) / I \ J-^T (C) \ rpT *(Cl \TT'^("^''^'^) / 



n^;^^^-^^(g, q', uj)=Y. ^.W(fi.)i^Jw*(^.')ni"r^(g, g'; c.), (46) 

JMW 

where Il{q,q']uj) is related with Yi{r,r']u) defined in Eq. (|^) as 

n:;(?^'^)(g,g» = rr2drr'2dr'j,(gr)n5?^'^)(r,r»>(gV). (47) 
Jo 



D. (vr + p + 5') model 



For the effective interaction ([T7|), the (vr + p + (7') model is commonly adopted. It gives 



with 



(2vr)= 



a,'--"^ ■ q 



a,^^') . q 



{2nf 



\^^(ab) . ^^(cd) 



(T&^'^ X q) ■ ((72^'"^^ X q 



In the grouped representation (pT|)- ([23|) , Wl and Wt are given [|^,|^,^ as 



mi 



gU^) + Tl{q,u;)Vl{q,u:) 



Lj2 _ g2 _ ^2 



/a//3 



9'a0i.q) + C'a/3r^^(9,^)r«(g,^) 



1^2 _ q,2 _ j^2 



(49) 
(50) 

(51) 
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where a and (3 denote or A, and and mp are the pion and the p- meson masses, 
respectively. The coefficient C^^ = -^^j-^^ is the ratio of the p-meson exchange coupling 



to the TT-meson one. The following vertex form factors are used, 



u^-q^-KV 



K{q.^) 



w2-g2_A2- 



(53) 



The Landau-Migdal parameters (?'(?) 's depend on q so weakly that we neglect their q de- 
pendence in the region of g < 3 fm~^ [^,^|. They are treated as free parameters. 
The effective interaction (|48|) is expressed in angular-momentum representation as 



Q 



where 

Wfir'\q,u^) = W,f,f\q,u) + lyrSr^lg,^), (55) 

with 

Wl^u''^ {q, uj) ^ anWL^'''''"'^ {q, u)aj,, (56) 
WTf,f\q,uj) = WT^''''''^''\q,uj)[Su^ - aj.aje]. (57) 

The coordinate representation of the effective interaction in Eq. (p^ is given by 



W\ff'^\ri,r2]uj) = - / q'^dq ji{qri)W^f/'^\q,u)ji'{qr2). (58) 

TT Jo 



oo 



The ring equation (p!8| ) reads as 

+ E rn'dnr^'dr^ n(°);;«^(r,ri;a;)iy^,^,^(ri,r2;a7)n:;,^,,(r2,r';^). (59) 

e. •'0 



E. Spin-longitudinal and -transverse response functions 

We define the spin-longitudinal and -transverse response functions as 

x<5[a;-(E„-Eo)], (60) 
Ip^'"''^ (g, ^) = ^ E ^ E ( ^0 I jTf ^) (g, u;) I vI/„ ) ( vl/„ I ^) t (g, ^) I VI/, ) 

x5[a;-(E„-Eo)]. (61) 
These are rewritten in the grouped representation as 

Rt'\q,uj) = -^-lmn,-^^'^^\q,uj), (62) 

A TT 

i?^t"^](g,^) = -iilmnT^["'31(g,c.), (63) 

A 71 

with the momentum diagonal parts of the polarization propagators ( |i5| ) and (^), 

UL'^"^\q,u;)^UL'^"^\q,q;u;) 

= 47r^(2J + l)aj,aj,,nif (g,g;u;), (64) 

Ur''^'^^\q;u)^W^T';^fKQ^q;u;) 

= 27r^(2J + - aMaje')Il'ju'\q,q;i^), (65) 

in 



where Eqs. (||), (||), (|38|) and (||) are used. 

As to the responses for nucleon probes, we assume that the ratio of the scattering am- 
phtudes, the A^A^ NA to the A^A^ — * A^A^, is f/^/fNi as is commonly done imphcitly |T3 
and exphcitly [^. Then relevant response functions are given by 



In \JnJ 
as was shown in Eq. (|I]). The uncorrelated ones are 



(66) 



(g, .;) = R'^r"\q, u) + ( ^ 1 Rt:^'^^ (g, u;) 



^{Q)u[NN\ . 
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(67) 



In the quasi-elastic region, the second term does not contribute since real A production does 
not occur. 

For (e, e') scattering, the cross section is expressed in the one-photon-exchange approxi- 
mation as 

.2\ 2 / 



dedfl 



0"A/ 



I SL{q,uj) + (tan 2^2^) ^^{q^oj) 



with the Mott cross section ctm, a transferred four-momentum q^ = {u, q), 



(68) 



and 



a scattering angle 9. The dynamic structure factor [|T^], Sl and St, are given byQ 

- (69) 

(70) 



SL{q, uj) = J2\{^n\ PcHq, ^) I ^0 ) 6[uJ - {En - ^0)], 
Sriq, ^) = E I ( I Jt^Q, uj) I ^0 ) f 5[uj - {En - E^)], 



where {pc, Jt) is the nuclear electromagnetic current operator. In this paper, we restrict 
ourselves to the transverse part Jt-, which is given by the sum of the one-body convection 
and magnetic currents, and the exchange current, Jt = J^^^ + Jt^^ + J'r'^^ ■ We neglect 



the convection current since its contribution is small Although the exchange current 
contributes to a certain extent, we do not take into account of this term since we only 
consider the responses of one-body operators. The magnetic current is given by 



rmag. 



2m N 



iqru 



{G'.'NNiqD + G!,Uiql) {^r') ,} (-.^^^^ x q) 



+ G,^^{ql){{rr'),{-^'''^'x<l 



(71) 



Some papers use the notation {Air /Mx)Sl^t instead of the present 5^^^, where stands for 
the target mass. 



1 1 



with the nucleon mass rriN. We take the magnetic form factors in the form of 



G.y a/sill) = /' 



7Q/3 



1 ^ 
A2 
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where A = 855 MeV/c |^ and fs are the magnetic strengths. Then we get 



2mN 

f-rNA r / (NA) 
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where 



f^NN — if^p + f^n) /2, /-JXat — if^p f^n) /2, 



7V 



(73) 



(74) 



with fip = 2.79 and Hn = —1.91. When we neglect the isospin-mixing and consider only 
T = target nuclei, the isoscalar and isovector parts do not interfere. Then we can neglect 
the isoscalar part since (/^^at/Z^jvat) ~ O-O^ is very small. Inserting Eq. (|73D into Eq. (|70D 
and using Eqs. ( |6TD and (|35|) , we get 

(75) 
(76) 

iq,u). (77) 
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where 



I^'''\q,u)=R'T=^'''''\q,uj) + 2 



u=0[NN] . 



jjU=0[N A] 

fIV 

J'yNN 



-JJV- 
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III. EFFECTS OF THE A-HOLE MIXING 

In the following two sections we present our numerical calculations of the isovector spin- 
longitudinal and -transverse response functions, i?^(g, a;) and Rrp{q,Lj), for the doubly (sub- 
)closed shell nuclei, ^°Ca, ^^O and ^^C, and analyze them from various points of view. 

In this section we compare Rj^x with and without A to see effects of A-mixing, and 
investigate relative importance of the process-decomposed response functions R^it- This 
manifests the A effects more clearly. We also discuss the energy-weighted and energy-non- 
weighted sums. 

Calculations are carried out by the ring approximation, the Tamm-Dancoff approxima- 
tion (TDA) and without any residual interactions. These results will be called the RPA, 
TDA and uncorrelated response functions, respectively. 
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A. Technical comments and choice of parameters 



Before presenting the numerical calculations, we make some technical comments and 
summarize the values of the parameters. 

In the previous section, H and \l/ include the center-of-mass motion, however, it is better 
to replace them by the intrinsic ones to isolate the structure part. Then the transferred 
energy u should be replaced by that to the intrinsic state, oj-mt — ^ — "-^recoil = — Q /{2AmN). 
Similarly the transferred momentum q is replaced by that to the relative motion between the 
active nucleon and the remaining {A — l)-nucleon core, gint = [{A — For simplicity 

we suppress the script "int" in Sects. [IIT| and [TV. 

We took the single-particle potential for N and A as 

U{r) = -iV + ^W) \ - V ' + ^^onl (78) 

1 + exp l^^j a p^i^ g^p (^)J 

with R = tqA^^^. VpQ^j is the Coulomb potential of the uniformly charged sphere with the 
radius parameter r^. The shape parameters are fixed to be tq = = 1.27 fm and a = 0.67 fm 
p^ . For the nucleon the spin-orbit potential depth Vis are fixed to be 6.5 MeV for ^^C, 
10.4 MeV for ^^O and 10.0 MeV for ^°Ca. The real potential depth V is so determined 
as to give the observed separation energy of the outermost occupied state. The imaginary 
potential depth W are fixed to be zero for the occupied (hole) states and 5.0 MeV for the 
particle states. For A we set V = 30 MeV and W = Vgs = 0.0 MeV since we do not have 
enough information for such virtual A appearing in the quasi-elastic region. 

The masses are chosen to be m^v = 940 MeV, = 1236 MeV, = 139 MeV, and 
rrip = 770 MeV unless explicitly mentioned. The coupling constants are fixed to be J^'^/Att = 
0.081, /a/Ztv = 2.00, and C^^ = 2.18. The cutoff parameters are set to be = 1300 MeV 



and Ap = 2000 MeV |13 



B. Effects of the A components 

Here we present the energy spectra of the response functions at q = 1.70 fm^^ for 
the charge exchange mode related with the (p, n)-like reactions {u = —1). From now on 
we suppress z/ on R. The Landau-Migdal parameters are taken to be {g'NNy dNAy 9'aa) = 
(0.6,0.4,0.5). 

In Fig. 1^, we show the RPA response functions R]^ of ^°Ca (a) without and (b) with A. 
The uncorrelated response functions Rl^^ are also shown, which are very close to each other 
for this LS closed shell. Slight difference comes from the spin-orbit force. Fluctuations seen 
at lower u are somewhat artificial because the spreading widths of the hole states are not 
included in the present calculation. 

15? 



When A is not included, Rj^ is enhanced below about 85MeV but quenched above that. 
On the other hand, Rj, is quenched below about 96MeV and enhanced above that. Conse- 
quently R]^ is softened but Rj. is hardened. 

Once A is introduced, Rj^ increases but Rj^ decreases for the whole quasi-elastic region 
as seen in Fig. |^(b). This is attributed to the coupling interaction W^^^\ which brings down 
the spin-longitudinal strength from the Ah to the ph sector but brings up the spin-transverse 
strength in the opposite direction. This is the essential effect of A. 

Fig. ^ shows the response functions with A for (a) ^^O and (b) ^^C. Qualitative features 
are common for all the nuclei and the RPA effects are stronger for larger A. The difference 
between R^^^ and R^^ is larger for ^^C than for other two nuclei, because ^^C is not the 
spin-saturated nucleus. 

To see the situation in more detail, we separate the response functions into the process- 
decomposed ones. The RPA response functions i?^ and Rrp of ^°Ca are decomposed into 
R^L^^ in Fig. 1(a) and R^^^^ in Fig. |(b). The main contribution comes from i?'^^', but 
that from R^^^'^ is also significant. The contribution from is negligibly small. An 

important point is that i?^^' and i?^^' contribute constructively, whereas i?^^^ and Rjf"^^ 
do destructively. This explains the shift of the spin-longitudinal and -transverse strengths 
seen above. 

To see the backward effect we show the TDA response functions without A in Fig. |^. 
Softening of i?^ and hardening of Rj^ are well developed, but their magnitudes do not change 
so much. Comparing with Fig. ^(a), we can say that the backward amplitudes in the ring 
approximation induce further enhancement of R^^ and quenching of Rrp. 

These behaviors are qualitatively understood in the following way. The formal solution 



of the RPA equation (29) is given by 

-1 



n 



i-nwi^l n(o), n(o) =ngi, + ngj,. (79) 



where ^^bk forward and the backward part of the uncorrelated polarization 

propagator. The fact to be kept in mind is that around the present momentum [q = 
1.7 fm~^) the spin-longitudinal effective interaction W^^\q,uj) is negative if g'^f^ < 0.7 but 
the spin-transverse one Wp''^\q,uj) is positive if g'^f^ > 0.25 for the present parameters (see 
Fig. 1). 

First let us consider the cases without A. For simplicity we treat 11 and W as c-numbers 
instead of matrices like in a Fermi gas model and consider a real single-particle potential. 
Then response functions are expressed as 

-2 



R 



n(°)iy ^R(°). (80) 



The real part of Uply changes the sign in the middle of the energy range concerned, negative 
at lower side of uj but positive at higher side, as is seen from the first term of r.h.s. of Eq. (pTI). 
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Hence HVI^PFt is negative but Hp^j^/iy^ positive at lower cj, and thus Rrp is quenched but 
i?^ enhanced if Wj^ is not so strong {Upl^Wj^ < 2). At higher u the situation is opposite. 
As the results the TDA response functions are softened for the spin-longitudinal mode but 
hardened for the spin-transverse mode. On the other hand H^bk is always negative (see 
Eq. (0)). Therefore R^^ is enhanced but Rj^ quenched by the backward amplitude in the 
whole energy region. 

Next let us consider the cases with A. In the first order of W the RPA response functions 
are given by 

j^lNA] ^ n(0)t^^V^,T[^^l4y^l. (81) 

The uncorrelated Ah polarization propagator n'^'^^'^^^ is real negative in the quasi-elastic 
region (see Eq. (pJ]) ). Therefore i?^^' is positive but is negative. This is the reason 

why i?^^' and i?^^^ contribute constructively, whereas Rj^^^ and do destructively. 



C. Sum rules 

Next we consider the energy-non-weighted and energy-weighted sums defined by 

^l,t(?) = / RLA<l^uj)du, (82) 
^l,t(?) = / ujR^^T{q,uj)duj. (83) 
When A is not included, the former behaves as 

Xl^iq)-.! (g-oo). (84) 

Here we used j = j for isovector currents since the isospin of the ground state is assumed 
to be zero. 

In the Fermi gas model there is the definite upper limit of the integral cUmax- In the case 
without A, it is 

2m m 

with the Fermi momentum pp- For the finite nucleus the upper limit extends to infinity 
in principle, because there is no sharp cutoff of the momentum distribution. However the 
response functions is small and decrease rapidly beyond c^max- Therefore we assume expo- 
nential damping beyond cu^ax to evaluate the integrations. 

In the case with A, R{q,uj) is sizable both in the ph and Ah sectors. For the sake of 
comparison, however, we make the same prescription as in the case without A, because we 
are interested in the strength distributed only in the quasi-elastic region in this paper. 



We present the energy-non-weighted sums X^y(g) in Fig. ^(a) for ^°Ca. We took pp = 
1.20 fm~^ to estimate cu^ax- The sums X° and X° of the uncorrelated responses are too 
close to distinguish. We found that they are also very close to the Fermi gas model value 



iQ[l-^) for Q<2 ^ ^^^^ 
1 for g>2 ' 



with Q = q/pF 1^^. For the uncorrelated responses reduction from unity reflects the Pauli 
blocking effect. 

The RPA correlation without A slightly increase X° but largely decreases Xj. as were 
seen in Fig. 0(a). Once A is included, is drastically enhanced, while Xj, is more quenched 
in the low q region but the quenching becomes smaller as q increases. 

In Fig. P(b), the energy- weighted sums X}^ j, are shown. We found that the RPA correla- 
tion without A hardly changes them from the uncorrelated cases, therefore we did not show 
the results for this case. Once A is introduced, the sum of the spin-longitudinal mode Xj^ is 
enhanced very much, whereas that of the spin-transverse mode X^ is only slightly affected. 

The sum rule says that the sums of the response function for an operator O are given by 

X\q) = {^o\dd^\^o)/A, (87) 
X\q) = {^o\[d,[H,d]]\^o)/A, (88) 

where O is assumed to be hermitian in Eq. (^). Pandharipande et al. calculated | \l'o ) 
exactly by using realistic nuclear force within the nucleon degree of freedom. Then they 
evaluated X° and X^ from these sum rules. It is found that their sum rule values are 
significantly larger than our results without A. Fig. |^(a) and (b) compare their values with 
ours for ^^O. For instance, their energy- non- weighted sums are about 19 % (longitudinal) 
and 10 % (transverse) and their energy- weighted sums are about 75 % (longitudinal) and 55 
% (transverse) larger than ours, at g = 1.70 fm~^. Such larger difference strongly indicates 
importance of the correlations beyond RPA. Note that in their calculation the effects of A 
are implicitly included in part through the A mediated three body force, but the processes 
expressed by B\^^^ and B\^^^ are not included. 



Thouless [^] proved that one gets the energy-weighted sum of the RPA response function 
by replacing I^E'o) in Eq. ( pSf ) by |\1'^^), the Hartree-Fock ground state wave function. 
This theorem explains why the RPA energy-weighted sums without A are very close to 
uncorrelated ones. It also supports that the difference between our results and those of 
Pandharipande et al. must be due to the nuclear correlations beyond RPA. This should also 
be reflected in the energy spectra of the response functions. The importance of the 2p-2h 
configuration mixing has been also pointed out by several authors [|19| , p2 . 



IV. DEPENDENCE ON EFFECTIVE INTERACTIONS 



In this section we investigate the effective interaction dependence of Rj^ rp and their 
process-decomposed components R^i^ ■ Some of the effective interactions are shown in Fig. ^ 
in the form of w4°t /('^)- Their [a/5] dependence comes only through g'^^- 

We present the g'j^^^ dependence of Rj^ and Rj- in Fig. H(a) and (b), and their g'^^ 
dependence in Fig. |^(c) and (d), respectively, for ^°Ca at g = 1.70 fm~^. We fixed g'j^^ = 0.4 
and g'^^ = 0.5 in the study of (^^^ dependence and g'j^j^ = 0.6 and g'^^ = 0.5 for g'j^^ 
dependence. The response functions so weakly depend on g'^^ that we do not discuss about 
it. 

Both Rp and Rp considerably depend on g'j^j^ as well as g'j^^. As g'j^j^ decreases, Rp 
becomes more enhanced at lower uj but less at higher u, while Rrp becomes less quenched 
at the lower side and less enhanced at the higher side. The g'j^^ dependence is more simply 
summarized. As g'j^^ decreases, both Rp and Rp increase, namely, Rp is more enhanced 
but Rrp less quenched. These features are qualitatively common for all nuclei ^''Ca, ^^O and 
^^C we analyzed. 

Such dependence is more clearly seen through the process-decomposed response func- 
tions. The g' dependence of R^i^ is shown in Fig. |10[ Since (?^^ controls the coupling 
strength between and A, R^^^^ is more sensitive to f?^^ than g'j^j^, whereas opposite is 
true for R^^^l 

Their behaviors are well understood by the interpretation given at the end of Sub- 
sec. [IIIB| . As g'j^j^ decreases the effective interaction Wl^^^'^ becomes more attractive as 
shown in Fig. |^. Consequently i?^^' is more enhanced at lower lj but more reduced at 
higher u. On the other hand, Wt^^^^ becomes less repulsive and therefore i?^^' is less 
quenched at lower u but less enhanced at higher u. These behaviors reflect in the g'^j^ 
dependence of Rp and Rrp seen in Fig. |^. 

As 

q'na decreases, Wl^^^^ becomes more negative but Wt^^^'^ does less positive. There- 
fore both i?f ^1 and 4^^^ are increased, consequently more enhancement of Rp and less 
quenching of Rp are resulted in. 

G.E. Brown and his collaborators |^7|^ advocated the scaling of effective masses of 
nucleon and mesons (except for pion) in the nucleus; e.g. m*p/mp ^ m*j^/mN. Correspond- 



ingly they claimed necessity of large g'j^j^ To see an implication of this proposal, we 



show in Fig. ^ the transverse response functions Rp with smaller rn*{= 0.75mp) and larger 
(?^^(= 0.8) together with the uncorrelated and the RPA responses with m* = nip and 
Q'nn = 0-6 as reference. We did not change the nucleon effective mass because our computer 
program cannot take into account its density dependence at the present and it must be the 
free nucleon mass at infinity. So the present calculation aims only to get feeling about the 
interaction dependence. 

Let us compare two cases, (a) a previous parameter set [m* = nip, {g'^ n i 9'n Ai 9'aa) 
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= (0.6,0.4,0.5)], and (b) the new one [m^ = 0.75m^, {g'^^^ 9n 9 aa) = (0.8,0.4,0.5)]. In 
Fig. ^(b) we also show Wt with smaller m*{= 0.75mp). It shows that Wt^^^^ is almost the 
same for the both cases around q ^ 1.70 fm~^ accidentally, and hence similar hardening 
is expected. However the positive Wt^^^'^ in the ) becomes very weak negative in 

the case (b), and hence i?^^' changes the sign. Consequently we see in Fig. |Tl| that the 
hardening of Rrp stays similar for the both cases but quenching is very much reduced in the 
latter. We must note that such change strongly depends on the momentum q as is seen 



in Figs, n and 15. We remark that Rj^ is also affected through the change of Its 
enhancement is reduced at lower u because of large g'j^^- 

Next we show the collectivity ratio Rj^/Rrp in Fig. [1^, for (a) ^°Ca and (b) ^^C at 
q = 1.70 fm^^. In the uncorrelated case, the ratio is, of course, close to unity, deviation 
from which is only due to the single-particle spin-orbit force. In the RPA calculation without 
A, the ratio is larger than unity at lower u but smaller at higher u. It is because R^ is 
enhanced at lower u and quenched at higher u, but Rrp behaves in the opposite way as was 
shown in Fig. 0(a). In the case with A, the RPA correlation drastically increases the ratio 
and makes it larger than unity almost for whole u. The case with {q'n ^ ^ 9'n a^ 9'aa) — (0-6) 
0.4, 0.5) gives larger ratio at lower uj than the almost universality case with (0.6, 0.6, 0.5) 
and changes more steeply as uj changes. Smaller g'jq^ brings down more R^^ strength from 
the A/i region. In the case with m* = 0.75mp, enhancement of R^ is suppressed due to 
larger g'^^^q = 0.8 , thus the ratio is smaller than the above two cases. 

The RPA effects became larger as the mass number increases. We note that these results 
of Rj^/Rrp cannot be compared with the (p, n) data at the present stage, because effects of 
distortion and absorption have not yet been considered. 



V. ELECTRON SCATTERING 



?(e,e') 



obtained by electron 



In this section we study the transverse response functions Rrp 
scattering (see Eq. ([77|)). 

To analyze the electron scattering, it has been known that we must take account of not 
only RPA correlation with A but also the mixing of 2p-2h or more complicate configurations 
p!9| , |22[] and the exchange currents []T9|,^, etc.. Here we do not intend to reproduce the 
experimental data, but want to see to what extent the RPA results depend on the magnetic 
transition ratio frfNA/ f^NN the effective interactions. 

As was mentioned in Sec. |, most of previous analyses used the ratio f-yNA/ f:^NN — 2-0. 
However, the SU(6) quark model ^0[ gives fryNA/ fyNN — 6-v^/5 ~ 1.70. Phenomenological 
analyses by Koch et al. and Kumano |Q gave the ratio 2.20 and 2.26, respectively. 
A reason for the discrepancy is that the ttN background scattering is renormalized in the 
phenomenological analyses but should be treated explicitly in the quark model |^ . 
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In Fig. |T3|, we compare the results with different values of f-yNA/ f^NN- The smaller 
the ratio the smaller the quenching, because the effect of A is essentially determined by 
the product ^{fjN/^/ f^^j^)Rlj^^\ We see some dependence on f-yNA/ J^nn when i?^^^ is 
sizable. 

We compare the results with various effective interactions for ^^C at q = 300 MeV/c 
in Fig. |n|(a) and at 400 MeV/c in (b), and for ^°Ca at 330 MeV/c in Fig. |15|(a) and at 
410 MeV/c in (b). Here we fixed f-yNA/f^NN = 2-20 Experimental dat a are shown as 
a reference. Compared with the energy spectrum of the uncorrelated case, the experimental 
spectrum is very much hardened, but the magnitudes are comparative. 

The RPA result with {q'n n ^ 9'n Ai Q'aa) = (0.6,0.6,0.5), which is practically the same as 
that of the universality ansatz (0.6, 0.6, 0.6), is quenched and hardened. As a result, the 
peak moves closer to the observed position, but the magnitude becomes much smaller than 
the data. 

The calculation with (0.6, 0.4, 0.5), which we used as the standard in Sec. |III|, places the 
peak at slightly higher energy and now at almost the right position. It also increases the 
magnitude though it is still smaller than the data. 

If we take m* = 0.75mp and g'^j^ = 0.8 as Brown and Rho suggested, the magnitude 
is increased very much and the peak comes very close to the experimental data. A good 



fit is seen in Fig. |T^(b) but overshooting in Fig. |15|(b). We must note that the good fit 
does not necessarily mean that the effective interaction is good because there must be other 
contributions. 

Qualitative features of our analysis are consistent with previous calculation of Alberico 
et al. , in which the RPA correlation gives reasonable hardening but underestimates the 
magnitude, the deficiency of which may be fulfilled by nuclear correlations beyond RPA (the 
2p-2h effects, etc.) and exchange currents. 

Since Rq^'^ ■* eminently depends on the effective interactions, it must be a good tool to 
discriminate them if the reliable estimation is possible of the other contributions such as 
2p-2h configuration mixing and exchange currents, etc.. Their estimation is beyond the 
scope of the present paper. 



VI. SUMMARY 

We studied the effects of the A-hole configurations on the spin-isospin response functions 
in finite nuclei in the quasi-elastic region. We removed the universality ansatz for the 
Landau-Migdal parameters and treated g'jyj^, g'jy^ and g'^^ independently. For this sake we 
formulated the response function method for a finite system consisting of and A in the 
ring approximation. 

We showed that the A-mixing is crucial for the enhancement of in the whole range 



1Q 



of the quasi-elastic region, and it promotes the quenching of If A is not included, i?^ 
and Rrp are both partially enhanced and partially quenched. We emphasize that reliable 
estimation of the effects of A is definitely needed for comparison with experimental data. 

Detailed analysis was carried out by dividing the response functions Rj^ j, into the process- 
decomposed ones, 4^1, R^? and R^? . The main contribution comes from i^'^l but 
contribution from -R^^' is also significant, whereas -R^^^' is negligible. The effects of A is 
mostly represented by R^it^- 

We showed that i?^^' and i?^"^' contribute constructively, whereas rI^^^ and i?^'^^ 
do destructively. This is the reflection that the negative interaction W^^^'^ between ph and 
Ah brings down the spin-longitudinal strength from the Ah to ph region, but the positive 
interaction wj,^^' brings up the spin-transverse strength in the opposite direction. Thus 
R^t^^ plays an important role for strong enhancement of Rj^, and i?^^' does some role for 
quenching of Rrp. 

Analysis of g' dependence of R^"'^'^ tells that R^^^'^ is very sensitive to g'j^/^ but not to 
9nn^ whereas is sensitive to g'j^j^ but not to g'j^fA- g'j^^ decreases, both R^ and Rj. 

increase, thus R^^ is more enhanced but Rj, less quenched. As g'jyj^ decreases, i?^ becomes 
more enhanced at lower u side but reduced at higher side, while Rj, becomes less quenched 
at lower side but more at the higher side. Consequently the choice of g'j^j^ and g'^^^ is 
crucial to determine the behavior of the response functions Rit, such as the collectivity 
ratio Ri^/Rrp. Thus we should relax the universality condition for g'^s. Effect of change of 
the p-meson effective mass is also presented. 

We further studied the interaction dependence of the isovector transverse response func- 
tions i?^'^ ■* for (e, e') scattering in comparison with the experimental data. The comparison 
must be useful to investigate the effective interaction at large q, though we need reliable 
estimation of the exchange currents, 2p-2h configuration mixing, etc.. 

In this paper we present the detailed analysis for ^"Ca. The similar analysis for ^^C is 



given in Ref. [54 



At the end we itemize some of the remaining problems we have to investigate. 1) For 
analysis of data of hadronic probes we must investigate reaction mechanisms such as distor- 
tions, multistep processes, etc.. This will be discussed in the forth-coming paper. 2) As we 
often mentioned nuclear correlations beyond RPA have to be evaluated. 3) We used Woods- 
Saxon shell model with the free nucleon mass. Hartree-Fock field should be used to keep 
consistency between the mean field and the residual interactions. 4) Nucleon effective mass 
with position dependence should be incorporated. 5) Spreading widths of ph propagation 
should be properly taken into account. At the present we only included that of the particles 
by an energy independent complex potential. 
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FIGURE CAPTIONS 



FIG. 1. Process-decomposed response functions. 

FIG. 2. Isovector spin-response functions, i?^ and Rj,, for ^'^Ca at g = 1.70 fm^^. (a) 
Without and (b) with A. The dotted and fuh hues denote and Rrp with RPA corre- 
lation, respectively. The dot-dashed and dashed lines represent R^^^ and r!^\ respectively. 
(.9nn^9na^9'aa) = (0.6,0.4,0.5) are used in RPA. 

FIG. 3. R^ and Rj, for (a) ^^O and (b) at g = 1.70 fm"^ with A. The notations of the 
lines and the values of g' are same as those in Fig. |2|. 

FIG. 4. Process-decomposed response functions, (a) R^^^^ and (b) i?^^', for ^^Ca at 
q = 1.70 fm^"*^. The dotted, dashed and dot-dashed lines denote R^^r^\ 2(/A//Ar)ii^^', and 
{f A / f n)"^ F^LT^ ^ respectively, appeared in the Eq. ([6^). i?^ and Rrp are shown by the full lines. 
The values of g' are same as those in Fig. ^. 

FIG. 5. Response functions for ^''Ca at g = 1.70 fm~^ without A calculated in TDA 
idNN ~ 0-6). The dotted and full lines show R^ and Rrp, respectively. The dot-dashed and 
dashed lines represent R^p^ and Rp \ respectively. 

FIG. 6. (a) Energy-non-weighted sums and (b) energy-weighted sums for ^'^Ca. The 
full line denotes the sums of the uncorrelated transverse response. The dotted line represents the 
RPA transverse results without A. The RPA transverse and longitudinal ones with A are shown 
by the dashed and dot-dashed lines, respectively. The uncorrelated longitudinal sum and the RPA 
longitudinal sum without A are very close to the uncorrelated transverse sum and not shown. 
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FIG. 7. (a) Energy-non-weighted sums and (b) energy-weighted sums for ^^O. The 
dashed and dot-dashed hues denote the transverse and longitudinal sums, respectively, calculated 
by Pandharipande et al. Our results without A are shown by the full (transverse) and dotted 
(longitudinal) lines. 

FIG. 8. Effective interactions Wr"^V(^^) at w = 80 MeV; (a) longitudinal and (b) transverse 
part. At g = fm^^ their values become equal to their g'^j^. The interactions with m* = rup and 
with m* = 0.75mp are denoted by the full and dotted lines, respectively. 

FIG. 9. The g' dependence of the response functions for '^'^Ca aX q = 1.70 fm~^. Left side: the 
g'j^j^ dependence of (a) Rj^ and (b) i?^ with g'j^^ = 0.4 and g'^^ = 0.5. g'j^j^ is set to be 0.5 (dotted 
line), 0.6 (dashed line) and 0.7 (dot-dashed line). Right side: the g'j^^ dependence of (c) and 
(d) Rrp with g'j^j^ = 0.6 and g'^^ = 0.5. g'j^^ is set to be 0.4 (dotted line), 0.5 (dashed line) and 
0.6 (dot-dashed line). The full line denotes the uncorrelated responses. 

FIG. 10. The g' dependence of the process-decomposed response functions for ^''Ca at 
q = 1.70 fm^^. The g'j^j^ dependence of (a) i?^^' and (b) Rj^^\ and the g'j^^ dependence of 
(c) R^^^^ and (d) rI^^^K The notations of lines and the values of g' are same as those in Fig. ^. 
The uncorrelated responses are not shown here. 

FIG. 11. Effective /o-meson mass dependence of the transverse response function for ^'^Ca 
at q = 1.70 fm^^ with the large (77V A^- '^^^ dot-dashed line represents the RPA result with 
m* = 0.75mp and (Sattv) SatA' 9aa) ~ (0.8,0.4,0.5). As reference, the uncorrelated response (full 
line) and the RPA one (dotted line) with (fi'jVAf' ^TVA' S'aa) ~ (0.6,0.4,0.5) are also shown. 

FIG. 12. Collectivity ratio R^/Rt for (a) ^°Ca and (b) at q = 1.70 fm"\ The thin full line 
denotes the ratio without the correlation. The dotted line represents the RPA results without A 
Wnn ~ ^■^]- '^^^ RPA results with A are shown by the dashed line [(^attv ^TVA' 5aa) ~ (0-6, 0.6, 
0.5)], the thick full line [{9'^^^ s'na^ b'aa) — (0-6) 0.4, 0.5)], and the dot-dashed line [m* = 0.75mp 
and (^W'SatA'Saa) = (0-8, 0.4, 0.5)], respectively. 



FIG. 13. Dependence on f^NA/fl^NN of St for ^^C{e, e') at q = 300 MeV/c. f^NA/f^NN is set 
to be 1.70 (dotted line), 2.00 (dashed line) and 2.20 (full line). {gNNidNAid'AA) is set to be (0.6, 
0.6, 0.5). 

FIG. 14. St for ""^^0(6,6') with various effective interactions at (a) q = 300 MeV/c and (b) 
400 MeV/c. fyNA/ f^fNN is set to be 2.20. The dotted and dashed lines denote the RPA results with 
{9'nni9'n Aid' aa) = (0-6, 0.4, 0.5) and (0.6, 0.6, 0.5), respectively. The dot-dashed line represents 
the RPA result with (ffjvAf > ffiVA' 9aa) = (0-^) O-^) O-^) — 0.75mp. The uncorrelated one is 

shown by the full line. The experimental data are taken from Ref. 

FIG. 15. St for ^''Ca(e,e') with various effective interactions at (a) q = 330 MeV/c and (b) 
410 MeV/c. The notations of the lines are same as those in Fig. |l^. The experimental data are 



taken from Ref. 115). 
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